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1 Introduction 

The presence of exotic quark matter among cosmic ray primaries has been a subject of 
experimental and theoretical studies for many years. It is safe to say that a milestone 
concerning this issue was the detection of the very intriguing “Centauro” events pQ, 
later interpreted as the explosion of a quark blob by Bjorken and Me Lerran [2], 
Over the years, many reports have been published and some retracted [3]. One of 
the latest examples of this class of events has been recently described by Basu et al. 
[3]. The importance of these phenomena for the knowledge of hadronic physics and 
its occurrence in nature can not be overstated. 

It is clear that the idea of quark blobs coming from outer space needs an appraisal 
in terms of production, ejection and survival of quark matter from its original site 
to the detectors. Since the last 30 years a stable form of that matter, the so-called 
strange quark matter (SQM) [5] has been studied and discussed. SQM, and its finite- 
A version (strangelets) is particularly suitable as a candidate for the primaries and 
expectations for its direct detection (for instance, in the AMS-02 experiment [6] have 
been published elsewhere [7j. However, there are some important caveats concerning 
the very process of fragmentation of SQM (assuming it is the true ground state of 
hadronic interactions), with the corresponding consequences for the total flux. We 
shall discuss below some of these issues, focusing on the physics of fragmentation. 
The goal would be to check whether a high enough strangelet flux is generated, inde¬ 
pendently of the acceleration processes [8j within the known models widely applied 
to nuclei. 


2 Strangelets from SS mergers, supernovae or both? 

Strange stars are the astrophysical realization of the SQM hypothesis. These should 
mimic ordinary neutron stars in many respects, and a lot of work has been performed 
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to distinguish them through some clear signature. This is far from easy, although for 
the moment we must assume that strange stars are present and form shortly after a 
collapse/supernova event |9j. 

The possibility that the very birth of strange stars ejects some SQM into the ISM 
has been considered long ago 0D33- The general picture relies on the development 
of instabilities ending in a fully turbulent propagation front HU, much in the same 
way as it is thought to happen in ordinary thermonuclear supernovae |T2]J. However, 
recent work [13j did not find any SQM ejection, at least within the explored physical 
conditions. Since (collapse) supernovae are very common, even a tiny fraction of bulk 
SQM would populate the galaxy with primaries and deserves consideration. 

In a paper addressing this question, namely what happens if bulk SQM is injected 
swiftly as a result of the propagation of a combustion (which may be entangled 
to cause of the supernova 0). we [H] have studied the fragmentation of SQM on 
their way out by collisions with oxygen nuclei using a simple spallation scheme [15] . 
Strangelets lose mass and kinetic energy in most cases, giving rise to a scaling law for 
the final mass mj? in terms of its initial value m; of the form 


( 'm F ^ _ 

mi 


/ A 0 E b \ 
' m p E 0i ' 






1 / 2 ( 23 — 1 ) 


( 1 ) 


where Ao is a scale for effective spallation of a fragment with that mass number, 
Eb is the binding energy released in the process, Eoi is the initial energy of an oxygen 
nucleus and D is a dimensionless quantity generally > 1 containing the radius of the 
strangelet and additional quantities. The bottom line of that calculation is that there 
is a real chance to obtain strangelets in the range A ~ 10 3 provided the spallation with 
oxygen nuclei is efficient enough. Of course, this does not address the injection mass 
spectrum, that is, the fragmentation prior to interactions with the oxygen nuclei, a 
subject that will be considered soon. 

An attempt to refine the crude model of Vucetich and Horvath [3] has been made 
in the last years [16]. Basically the model replaced the empirical Boyd-Saito [15] with 
an abrasion-ablation scheme of the type usually employed in the GEANT4 algorithm 
and similar procedures. We found that at around E/A ~ 10 MeV/A the spallation 
(abrasion) quickly dominates the fusion and scattering cross-sections, giving rise to 
a much moderate downsizing of simulated strangelet-nitrogen collisions. Therefore, 
it is likely that the first models overestimated the reprocessing of strangelets and 
A ~ 10 2 — 10 3 fragments could be produced only if the initial baryon number was not 
very high. Thus, the initial mass distribution becomes an essential ingredient. 

The idea of fragmentation of bulk ejected SQM itself was never considered in 
depth because the binding energy of the former always increases with baryon number 
A. At first sight, producing fragments out of a big chunk would go against energetics 
because the process must work “uphill”. Madsen H7j equated the gravitational energy 
to the surface energy to estimate that around A ~ 10 38 the SQM should break up in 
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chunks of that A. His picture of big nuggets orbiting a SS lead him to identify orbital 
collisions as the mechanism to further fragment the chunks into A ~ 10 3 strangelets, 
that is, a very substantial reprocessing indeed. 

In complete analogy with bulk nuclear matter, it is expected, however, that a 
consideration of the free energy of the system (not just the binding energy) indicates 
the fragmentation point and also the mass distribution of the fragments. The physi¬ 
cal situation of any ejected SQM in bulk corresponds to an expansion and cooling of 
the system, and therefore excluded volume at the fragmentation point and vacuum 
“melting” should be important as well. Note that this is independent of reprocessing 
by collisions. Moreover, such a situation would be encountered in the first moments 
after a supernova shock breakout, but also in the situation of ejection at the merger 
of two compact strange stars [18]. The latter process could eject as much as 1O _4 M 0 
in the form of SQM, but this would produce either planetary mass chunk(s) or re¬ 
quire very efficient reprocessing, as explained above, to later match the reported CR 
primary range. We shall now present the results of a calculation intended to explain 
the fragmentation and the mass distribution of the fragments, without considering 
further collisions or acceleration. 


3 Multifragmentation of SQM 


Statistical multifragmentation models are known to reproduce the fragment distri¬ 
bution of ordinary heavy nuclei and can be applied whenever the temperature T is 
of the order of the binding energy. In the grand canonical ensamble, the partition 
function reads [T9] 

= v (—^ e - (2) 
with f A the internal free energy of the fragment 

f A = -W A + oA 2/3 + CA 1/3 , (3) 

representing the bulk, surface and curvature terms respectively. Performing the 
Fourier transform of eq.(2) Z, and recalling that 
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the problem of the fragmentation is reduced to finding the singularities of the 
isobaric partition function, yielding the gas and liquid pressures as 
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where v = fi + W is the (shifted) chemical potential. The fragmentation spectrum is 
obtained from the derivatives of the pressure assuming chemical equilibrium. How¬ 
ever, when these quantities were calculated for SQM within the MIT bag model 
framework and in the color-flavor-locked (CFL) state [20], the surprising result that 
the peak of the fragments happens in the region A <Cl was obtained in all cases, even 
introducing a simple Van der Waals approximation. This may be due to the hidden 
assumption about a constant energy of the vacuum (represented by an MIT value B. 
However, physical considerations can be made in favor of a “melting” of the vacuum 
picture: a fraction of the vacuum energy is being used to create the surfaces and 
curvature energies, therefore there should be an energy density difference between 
the phases 


fh) — W'o + Bvu q , 

Wg = Wq + BVgaS- 

After this correction, the distribution shifts to higher values of A. We have normal¬ 
ized the probabilities to a total of 1O 4 M 0 as expected from the above considerations 
M- The masses of the fragments (strangelets) reach higher values when progressively 
higher temperatures are considered (Fig. [[]). However, as is well-known the stability 
of the strangelets diminishes with rising temperature, and an increasing number of 
fragments should decay into ordinary hadrons 

4 Conclusions 

From the calculations performed in a series of studies we may conclude that either 
SQM does not like to fragment, or if it does an overwhelmingly large fraction of the 
original baryon number would be ultimately go into ordinary particles, not strangelets. 
A similar calculation by Biswas et al [22] has pointed out the sensitivity to the treat¬ 
ment of the binding energy dependence, while the results above become odd without 
a substantial role of the vacuum melting. If the fragmentation destroys most of 
strangelets, the flux of exotics in CRs could be minuscule. Actually, the (optimistic) 
expression for the latter given by Madsen [23] 

F — 5 x 10 5 (m 2 yr sterad)- 1 x AL 4 x M_ 2 x V m \ x t 7 , (7) 

with R _ 4 the rate of mergers in units of 10~ 4 yr~ 1 , M_ 2 the ejected mass going 
into strangelets in units of 1O~ 2 M 0 should be multiplied by an efficiency factor of at 
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Figure 1: The distribution function of the fragments (strangelets) normalized to 
10~ 4 M o for three different values of the temperature (indicated). Fragmentation 
at higher temperatures favor a tail of heavier particles, but also destabilizes the small 
A strangelets. As a result, the net number of surviving strangelets decreases at higher 
temperatures. 

least 10 -5 representing the strangelet fragility under decays. This is independent of 
the possibility that the ejection could indeed be zero (pSj) in these events. Adding 
the caveats already made in the supernova case, we conclude that if strangelets are 
positively detected after all, the fragmentation should proceed out of equilibrium, a 
possibility yet to be explored. 

Acknowledgements 

JEH expresses his thanks to the organizers of the CSQCD IV conference for the real¬ 
ization of a live and productive event in which many ideas and paths were discussed 
and suggested to him. The financial support of the Fapesp (Sao Paulo State) and 
CNPq (Brazil) agencies is gratefully acknowledge by both authors. 


References 

[1] C. Lattes, Y. Fujimoto and S. Hasegawa, Phys. Rept. 65 (1980) 151. 

[2] J. Bjorken and L. Me Lerran, Phys. Rev. D 20 (1979) 2353. 

[3] M. Ichimura et ah, Nuovo Cimento A 106 (1993) 843. 


5 







[4] B. Basu et al., Astropart. Phys. 61 (2015) 88. 

[5] A. Bodmer, Phys. Rev. D 4 (1971), 1601; H. Terazawa, Tokyo U. Report INS-336 
(1979); E. Witten, Phys. Rev. D 30 (1984) 272 

[6] See the homepage http://www.ams02.org/ for a complete report n the AMS-02 
experiment. 

[7] J. Madsen, Phys. Rev. D 71 (2005) 014026. 

[8] G.A. Medina Tanco and J.E. Horvath, Astrophys. J. 464 (1996) 354. 

[9] O.G. Benvenuto and J.E. Horvath, Phys. Rev. Lett. 63 (1989) 716. 

[10] O.G. Benvenuto and J.E. Horvath, Mod. Phys. Lett. A 4 (1989) 1085. 

[11] J.E. Horvath, Int. Jour. Mod. Phys. D 19 (2010) 523. 

[12] C.M. Malone et ah, Astrophys. J. 782 (2014) 11. 

[13] M. Herzog and F.K. Ropke, Phys. Rev. D 84 (2011) 083002. 

[14] H. Vucetich and J.E. Horvath, Phys. Rev. D 57 (1998) 5959. 

[15] R.N. Boyd and T. Saito, Phys. Lett. B 298 (2003) 6. 

[16] L. Paulucci and J.E. Horvath, Phys. Lett. B 733 (2014) 164. 

[17] J. Madsen, Lecture Notes in Physics 516 (1999) 162. 

[18] A. Bauswein et al., Phys. Rev. Lett. 103 (2009) 011101. 

[19] K.A. Bugaev et al. Phys. Rev. C 62 (2000) 44320. 

[20] M. Alford et al. Phys. Rev. D 64 (2001) 074017. 

[21] L. Paulucci and J.E. Horvath, Jour. Phys. G 36 (2009) 095202. 

[22] S. Biswas, et ah, Phys. Lett. B 715 (2012) 30. 

[23] J. Madsen, Jour. Phys. G 28 (2002) 1737. 


6 



